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ABSTRACT 

We show that the hierarchical assembly of cold dark matter (CDM) haloes preserves 
the memory of the initial conditions. Using N-body cosmological simulations, we 
demonstrate that the present-day spatial distribution and kinematics of objects that 
formed within early (z > 10) protogalactic systems (old stars, satellite galaxies, globu- 
lar clusters, massive black holes, etc.) depends mostly on the rarity of the peak of the 
primordial density field which they originally belonged to. Only for objects forming 
at lower redshifts the exact formation site within the progenitor halo (e.g. whether 
near the center or in an extended disk) becomes important. In present-day haloes, 
material from the rarer early peaks is more centrally concentrated and falls off more 
steeply with radius compared to the overall mass distribution, it has a lower velocity 
dispersion, moves on more radial orbits, and has a more elongated shape. Population 
II stars that formed within protogalactic haloes collapsing from ^ 2.5a fluctuations 
would follow today a r~ 3 - 5 density profile with a half-light radius of 17 kpc and a 
velocity anisotropy that increases from isotropic in the inner regions to nearly radial 
at the halo edge. This agrees well with the radial velocity dispersion profile of Galaxy 
halo stars from Battaglia et al. (2005) and with the anisotropic orbits of nearby halo 
stars. 

Key words: methods: N-body simulations - Galaxy: halo - kinematics and dynamics 
- galaxies: formation - haloes - star clusters 



1 INTRODUCTION 

In a Universe where cold dark matter (CDM) dominates 
structure formation, the haloes of galaxies and clusters are 
assembled via the hier archical merging an d accretion of 
smaller progenitors (e.g. lLacev fc Colelll99ot) . This process 
causes structures to relax violently to a new equilibrium 
by redistributing energy among the collisionless mass com- 
ponents. Early stars formed in these progenitors behave as 
a collisionless system just like the dark matter particles in 
their host haloes, and they undergo the same dynamical pro- 
cesses during subsequent mergers and the buildup of larger 
systems like massive galaxies or clusters. It is of crucial im- 
portance in galaxy formation studies to explore the efficiency 
of the mixing process and see if any spatial or kinematical 
signatures exist in material that collapses at different epochs 
and within peaks of the primordial Gaussian density field of 
different rarity. 
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In this paper, we use a suite of high-resolution cosmo- 
logical N-body simulations to analize the distribution and 
kinematics within present-day galaxy haloes of dark mat- 
ter particles that originally belonged to selected branches 
of the merger tree. These properties are particularly rel- 
evant for the baryonic tracers of early CDM structures, 
e.g. the old stellar halo which may have originated from 
the disr uption at high reds hift of numerous dwarf proto- 
galaxies feullock et al.l200Cfl. the old halo globular clusters, 
and also giant ellipticals Jcao et al.ll2004) . The end prod- 
uct of the entire merger tree is a triaxial cuspy dark mat- 
ter halo (|Dubinski fc Carlberd Il99lt [Navarro et alJ Il996t 
iMoore et al .1 119991: iDiemand et alJl20uah a small fraction of 
early progenitor systems survive the merging process an d 
end up as dark matter substructures llGhigna et aljlT99ct) . 
Since rare, early ha loes are strongly bia s ed towards over- 
dense regions (e.g. ICole fc Kaiserf jl989l ISheth fc Tormenl 
1999), i.e. towards the centers of larger scale fluctuations 
that have not collapsed yet, we might expect that mate- 
rial originating from the earliest branches of the merger 
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tree is today much more centrally concentrated than the 
overall halo. Indeed, a "non-linear" pea k biasing has been 
discussed previously by several auth ors ijMoore et, al]ll99ji 
IWhite fc Springell200a lMoorell200lh . 

Here we show that the distribution and kinematics of 
"old material" within present-day galaxy haloes depends pri- 
marily on the rareness of the peaks of the primordial den- 
sity fluctuation field it originally belonged to. Specifically, 
today's properties of objects that formed in old rare den- 
sity peaks above vaiM, z) [where cr(M, z) is the rms fluctu- 
ation in the density field linearly extrapolated to redshift 
z smoothed with a top-hat filter of mass M. A "one a 
peak" corresponds to the characteristic mass "M*(z)" and 
cr(M«, z) ~ 1.69], depend largely on v and not on the partic- 
ular values of z and M. Such centrally concentrated compo- 
nents are isotropic in the inner part, just like the host galaxy 
haloes, but rapidly become more radially anisotropic further 
out. The plan of the paper is as follows. In § 2 we describe 
the numerical simulations and how to define and trace high- 
er particle subsets. In § 3 we analyse the present-day distri- 
bution of these subsets. We derive a simple empirical fitting 
formula for p(r, v\ the mass density profile of all progenitors 
above va, which approximates the results of our N-body cos- 
mological simulations for 1 < v < 4. § 4 discusses the impli- 
cations of our findings for old stellar populations. We argue 
that such centrally concentrated components are predicted 
to be isotropic in the inner part, just like the host galaxy 
halo, but to rapidly become more radially anisotropic further 
out. This is quantitative with t he radial velocity dispe rsion 
profile of Galaxy halo stars from lBattaelia et alJ l|2005r) and 
with the anisotropic orbits of nearby halo stars (/3 ~ 0.5, 
IChiba fc BeerslEoOOl) . Finally, we present our conclusions in 
§ 5. 



2 METHOD 

We identify collapsed high-er peaks at different epochs 
within high-resolution cosmological N-body simulations, 
mark them, and analyse the distribution and kinematics of 
this material at redshift zero. Details about the simulations 
and halo finding method are given in the next two subsec- 
tions. 

2.1 Simulations 

The simulations have been performed using PKDGRAV, 
a parallel N- body treecode written by Stadel and Quinn 
(Stadcl 200lT): c osmological and numeri cal parameters are 
the same as in iDiemand etafl J2004cft . The present-day 
haloes that we analyse are labeled "D12" (cluster) and "GO" 
to "G3" (galaxies). We also analyse an additional smaller 
galaxy halo (M vir = 10 11 M ), labeled "G4". 

Figures 1 and 2 show the density field of the region con- 
taining the G0-G4 galaxy haloes. Most of the material be- 
longing to high-a peaks in the high-resolution region can be 
found in the four most massive galaxy haloes, while smaller 
structures tend not to have any progenitors which meet the 
selection criteria. G4 is a special, early-forming small galaxy 
halo (with a high concentration of c = 18) that contains sim- 
ilar fractions of high-a material as the more massive galax- 
ies. From Figure^it is clear that G4 is not a representative 



Table 1. Present-day properties of the six simulated dark mat- 
ter haloes. The columns give halo name, spline softening length, 
number of particles within the virial radius, virial mass, virial ra- 
dius, peak circular velocity, and radius to the peak of the circular 
velocity curve. 



Halo 
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M vir 


' vir 


K,max 


r V c ma - 




[kpc] 


[10 6 ] 


[10 12 M Q ] 


[kpc] 


[km/s] 


[kpc] 


DV2 


1.8 


14.0 


305 


1743 


958 


645 


GO 


0.27 


1.7 


1.01 


260 


160 


52.2 


G*l 


0.27 


1.9 


1.12 


268 


162 


51.3 


G2 


0.27 


3.8 


2.21 


337 


190 


94.5 


G3 


0.27 


2.6 


1.54 


299 


180 


45.1 


G'l 


0.27 


0.25 


0.144 


138 


96.4 


15.0 



system for its mass range. All "G" haloes form in the one 
region of ~ (10 Mpc) 3 that was selected from a (90 Mpc) 
box and re-simulated at high resolution. The region has a 
mean density of 0.5 times the background density at z — 0. 
The properties of the six haloes are given in Table 1. 

The cosmological parameters are (fi m , Qa, as, h) = 
(0.268,0.732, 7,0.71) . The value of a 8 = 0.9 given in 
IDiemand et alJ <l2004cf) is not correct: we found that, due to 
a mistake in the normalization, our initial conditions have 
less power than intended. According to linear theory, mass 
fluctuations grow proportional to the scale factor in a flat 
Q m = 1 Universe, which is a good approximation to the 
adopted ACDM Universe at z > 2. Therefore the scale fac- 
tor at collapse of a given halo in our simulations would be 
0.78 times smaller in a as = 0.9 model. Throughout this pa- 
per, together with the collapse redshifts given for our low-ers 
simulations, we will also state the rescaled collapse redshifts 
in a as = 0.9 Universe using the exact growth function for 
a (f2 m ,fiA) = (0.3,0.7) cosmology (see table |2Jl. Using a 
as = 0.9 resimulation of the "D" cluster we confirm that 
the 2 = results (in units of the scale radius) are consistent 
with those of the <ts = 0.7 simulation when we use the Z0.9 
outputs of the as = 0.9 run to select progenitor haloes. 

2.2 Tracing progenitor material 

Since baryonic objects (i.e. early protogalaxies, first stars) 
will form in the inner parts of their host gravitational po- 
tential, it is important to assess how different their z — 
distribution is when only the inner part of a progenitor halo 
is traced. We will therefore compare the final distribution of 
material marked within the central regions of early haloes 
with the entire marked progenitor. We will show in this sec- 
tion that there is no difference (except of course in the traced 
mass fraction) if the traced objects are early (z > 10) pro- 
togalaxies. 

Progenitor haloes a re identified using the friends-of- 
friend algorithm (FOF) jDavis et alJ 119851) . For simplicity 
we use a fixed linking length b — 0.164 at all output times. 
Our results are not sensitive to this choice. For example 
at redshift 10.5 we found 52 haloes more massive than 
4.35 x 10 9 Mq in the high resolution region of the galaxy 
simulation. When we use a linking length two or three times 
smaller and mark the linked material in the same 52 groups 
we obviously find a smaller fraction of marked matter but 
very similar density profile shapes and kinematics (see Fig. 
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Figure 1. Density map of the high-resolution region of run G at z = 13.7. FOF groups more massive than 4.9 X 10 7 Mq (=84 mrjM) 
are marked in green color. 
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Figure 2. Density map of the high-resolution region of run G at z = 0. The marked particles (in green) were selected at z = 13.7 in 
groups above 4.9 X 10 7 M , see Fig. [I] The squares enclose the virial radii of the five galaxy haloes analysed (GO to G4). 
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Figure 3. The present-day distribution and kinematics of dark matter particles selected at z = 10.5, averaged over four Milky Way-size 
simulated haloes. We marked the 52 most massive groups found with FOF using different linking lengths: b = 0.164 (dotted lines), 
b = 0.164/2 (long-dashed lines) and b = 0.164/3 (short-dashed lines). Particles selected at the same redshift with overdensities above 
1000 have similar present-day distributions (dash-dotted lines). For comparison we also plot the same quantities using all particles from 
the final structure (solid curves). 



. We also selected the cores of haloes at the same redshift 
by marking all particles which have a local density (cal- 
culated from an SPH kernel over 32 nearest neighbours) 
10 3 times higher than the mean matter density and ob- 
tained similar z — distributions. This demonstrates that at 
redshift zero, particles originating from the cores of high-er 
haloes are distributed in the same way as all the progenitor 
halo material, and justifies our choice of using all halo parti- 
cles as tracers for baryonic objects which would likely form 
in the very center of their hosts. A large number of tracers 
give two important advantages. First, it allows us to reliably 
estimate density, velocity dispersion, and anisotropy profiles. 
Second, it makes the results more robust against numerical 
effects: since two-body relaxation completely changes the 
orbits of many individua l particles, even in hig h-resolution 
cosmological simulations (iDiemand et al . 2004a), using only 
one most bound particle as a tracer for (say) a Population 
III remnant is not a safe choice. 

We have shown above that the final distribution and 



kinematics of particles from high-redshift, low-mass progen- 
itor haloes is insensitive to their original location within the 
host potential. The orbital energy associated with the merg- 
ing of these small subunits into larger systems determines 
the present-day distribution and dominates the relatively 
small differences in the orbital and potential energies of par- 
ticles within their hosts. This is not true at lower redshifts 
and for progenitor masses closer to the mass of the z — 
parent halo. For example, progenitors above lo~ marked at 
z = 1.6 (i.e. with masses above 7.4 x 10 9 M©) show only 
mild radial and velocity bias when the entire FOF groups 
are marked (Fig. ■ the difference from the total dark mat- 
ter distribution becomes much larger, however, if we select 
within these haloes only particles with overdensity greater 
than 10 4 , as this material is more concentrated, colder, and 
on more radial orbits (just like in the low-mass 2a selection, 
see Fig. 0. 
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Table 2. Properties of 4, 3.5, 3, 2 and lcr fluctuations collapsing at different output rcdshifts. Here Ng IonpB 
is the number of groups found with mass > M, and / is the fraction of the overall mass in the high-resolution 
region which is inside these groups (the total mass in the high-resolution region is 3.6 X 10 13 Mq in run G, 
and 1.6 X 10 15 M in run D12). The fractions differ from the average value, / = er(c(v/\/2), as the high- 
resolution regions analysed are relatively small and are not representative patches of the Universe. We did 
not use outputs when the mass scale lies below ten times the particle mass, or when the number of groups 
is smaller than two per present-day halo. In the cluster run (D12) we have only analysed one representative 
case for lcr and 2a progenitors. 



V 


z 


^0.9 




-^groups £ G 


fee 


^groups £ .D12 


/eD12 








[Mq] 










4 


10.5 


13.8 


3.4 x 10 9 


- 


- 


24 


1.2 x 10~ 4 


4 


8.7 


11.5 


1.8 x 10 10 


- 


- 


14 


2.9 x 10~ 4 


4 


7.4 


9.8 


6.2 x 10 10 


- 


- 


8 


5.8 x 10~ 4 


4 


5.7 


7.7 


3.6 x 10 11 






3 


1.2 x 10 — 3 


3.5 


16.3 


21.2 


5.8 x 10 6 


72 


2.5 x 10~ 5 


- 


- 


3.5 


13.7 


17.9 


4.9 x 10 7 


35 


1.3 x 10" 4 


- 


- 


3.5 


10.5 


13.8 


8.4 x 10 8 


(5) 




236 


3.4 x 10~ 4 


3.5 


8.7 


11.5 


5.2 x 10 9 






126 


1.0 x 10~ 3 


3.5 


7.4 


9.8 


1.9 x 10 10 






60 


1.7 x 10~ 3 


3.5 


4.4 


6.0 


6.6 x 10 11 






8 


8.5 x 10" 3 


3 


13.7 


17.9 


5.9 x 10 6 


670 


3.0 x 10" 4 






3 


10.5 


13.8 


2.2 x 10 8 


52 


8.7 x 10" 4 


1659 


7.1 x 10~ 4 


3 


8.7 


11.5 


1.0 x 10 9 


23 


1.6 x 10~ 3 


1205 


2.5 x 10~ 3 


3 


7.4 


9.8 


4.3 x 10 9 


8 


2.0 x 10~ 3 


501 


4.3 x 10~ 3 


3 


4.4 


6.0 


2.3 x 10 11 






27 


1.4 x 10~ 2 


2.5 


10.5 


13.8 


1.8 x 10 7 


1426 


2.5 x 10~ 3 






2.5 


8.7 


11.5 


1.0 x 10 8 


558 


5.1 x 10" 3 






2 


8.7 


11.5 


7.6 x 10 6 


12740 


1.1 x 10~ 2 






2 


7.4 


9.8 


4.8 x 10 7 


2657 


1.5 x 10~ 2 






2 


4.4 


6.0 


6.3 x 10 9 


43 


2.1 x 10"2 


2459 


2.1 x 10~ 2 


1 


2.4 


3.5 


4.8 x 10 s 


1413 


0.14 






1 


1.6 


2.4 


7.4 x 10 9 


108 


0.13 


4258 


0.24 


1 


1.1 


1.8 


3.9 x 10 1(1 


24 


0.11 






1 


0.8 


1.5 


1.0 x 10 11 


9 


0.11 







3 DISTRIBUTION AND KINEMATICS OF 
HIGH-cr PEAKS 

We use the linear growth approximation to structure for- 
mation to calculate the masses of 4, 3.5, 3, 2 and lcr fluc- 
tuations collapsing at a given redshift. The redshifts of the 
simulation outputs and the halo masses corresponding to 
these fluctuations are given in Table H These values are for 
the cosmological model we have simulated (O m , Oa, <js, h) = 
(0.268,0.732,0.7,0.71). In a a 8 = 0.9 the same fluctuations 
would collapse earlier (at Zo.g, see S I2.lt . 



3.1 Radial distribution 

High-cr material is strongly biased toward the center of 
present-day haloes (see upper panels of Figs. l4lto ll2L with 
the rarer peaks showing stronger bias. Figures 2] to [5] show 
different (M,z) selections corresponding to the same value 
of v. While the shapes of these profiles are similar for a 
given v, the normalisation (or traced mass fraction) gener- 
ally grows with increasing mass threshold and to lower red- 
shifts. The density profiles are averaged over the four galaxy 
haloes of similar mass. The scatter from halo to halo in their 
total density profiles and that of their subsets is relatively 



small. But there are substantial halo-to-halo variations in 
the anisotropy parameter /3(r), which we discuss in § 13.21 

Mass densities are fitted with a general Q/37-profile that 
asymptotes to a central cusp p(r) oc r~ 7 : 



p{r,v) = 



(1) 



(r/n,)T[l + (r/rv^p-T)/" 

For comparison the NFW jNavarro. Frenk fc Whitel | l996l) 
profile has (a, (3, 7) = (1, 3, 11. while the lMoore et1flT<ll999h 
profile has (a,f3, 7) = (1.5,3,1.5). We fix a = 1 and the 
inner slope to 7 = 1.2, which is the best-fit slope for the 
7312 cluster w hen resolved at very high resoluti on (ttidm = 
3.0 x 10 5 M Q ) JPiemand. Moore fc Stadell2005h . We fit the 
entire dark halo using an outer slope of /3 = 3 to determine 
the scale radius r s = r V i r /c, where c is the concentration. To 
approximate the high-cr subset profiles we use a smaller scale 
radius r v = r s / f u (corresponding to a higher concentration 
c„ = /„c) and also a steeper outer slope f3 v . The /„ and /3 V 
values used in the plots are calculated with simple empirical 
formulae which approximately parameterise the entire range 
of profiles, i.e. peaks above 1 to 4a and haloes ranging from 
a low concentration (c = 4.5) cluster halo to a small, c = 17 
galaxy halo: 

.1.6 



,/f* , U = exp^/2) , /3 V = 3 + 0.26V 1 



(2) 
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Figure 4. The present-day distribution and kinematics of peaks above la averaged over four parent Milky Way-size galaxies. Here and 
in subsequent plots, the same-cr peaks selected at the highest and up to the next three output rcdshifts (given in Table 2) are plotted as 
dash-dotted, dotted, long-dashed and short-dashed lines, respectively. For comparison we also plot the same quantities using all particles 
(solid curves). The triangles in the upper left and lower right panels show the empirical fitting functions in J2J and IsV . 



The values for the 1 to 4a peaks are given in Table 01 and 
the profiles are plotted in the upper left panels of Figures 
H to E2 with open triangles. The fits are just approximate 
but they are reproduce well at least one profile for each of 
the v values. Other parameters and functional forms could 
fit some of the data better. 

For easier comparison with extragalactic observational 
data of old stellar or globular cluster populations we also fit 
our high-o subset profiles with a deprojecte d R 1 ^ law, us- 
ing th e accurate numerical approximation of lMarauez et all 
(2001): 



p(r) 

P 
b 



-Pg-Mr/fle) 1 / 4 



po{r/Re 

1.0 - 0.6097(1/4) + 0.05563(l/4) 2 
1.9992 x 4- 0.3271 . 



(3) 



The best-fit effective radii R e are given in Table El Like the 
TV values they scale with r s and not with the virial radius. 



3.2 Kinematics 

The lower left panels of Figures 21 to 1121 show the three di- 
mensional velocity dispersion profiles of D = a, 2 + of , where 
o r and (j t are the radial and the tangential velocity dis- 
persions (of = a 2 + o" 2 ,). The lower right panels show the 
anisotropy parameter /3 = 1 — 0.5o 2 /o 2 as a function of 
radius. The high-o- subsets are clearly slower and on more 
radial orbits than the entire dark matter component at the 
same radius. As with the density profiles, the differences de- 
pend on v, i.e. the highest o peaks have the lowest 030 and 
the largest (3 values. 

Like the density profiles also the velocity dispersions 
are averaged over the four galaxy-size parent haloes. There 
is little halo-to-halo scatter in p(r) and ow{r), but there is 
substantial scatter in the anisotropy parameter (3(r), espe- 
cially near the virial radius. We have checked that the radial 
bias exists in each data set and is not just on average. The 
variation in f3 over the four galaxies is about 0.1 within 10% 
of the virial radius, for both the parent halo and the progen- 
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Figure 5. Distribution and kinematics of peaks above 2a in present-day galaxy-size haloes. Line styles are as in Fig. [I] The squares in 
the lower left panel are calculated from the Jeans equation J1J using approximations for p(r) and /3(r) as input (plotted with squares in 
the upper left and lower right panels). Dispersions from an isotropic model (filled circles) and a model with constant radial anisotropy 
(3(r) = 0.5 (open circles) do not fit the corresponding data (long-dashed lines) very well. 



itor subsets. Further out the scatter becomes much larger: 
particles in the parent halo (subsets) have a total spread of 
0.45 (0.35) near r vir . 

3.2.1 Jeans equations 

The origin of the smaller velocity dispersions of the subsets 
can be traced to the fact that high-cr peaks form closer to the 
main progenitor of the present-day parent halo, hence they 
join the system with small infall velocities and at early times 
when its potential well is much shallower. Since they have 
been part of the parent halo for a long time, they are likely 
in dynamical equilibrium with the host, i.e. their density and 
velocity profiles should be a stationary solution to the Jeans 
equation in the potential, ^(r), of the parent halo. Setting 
time derivatives to zero and assuming sph erical symmetry, 
the J eans equation can be written as (Bi nnev fc Tremaind 
1987) 



Tr ^ 



2/3 a d* 
H po r + p—— — 

r ar 



(4) 



Approximating the anisotropy 
iDiemand et"al] (l2004d) . 

/3(r)=/Mr/r vir ) 1/3 , 

equation (QJ admits the solution 



parameter as in 



(5) 



p(r)a r (r) = Ge 



-6/3(r) 



y p(y)M{y) dy 



(6) 



Figure |S] depicts this solution for 2a material. We use equa- 
tion |IJ with (a,f3„,j) = (1,5.4, 1.2) and r v = r a = rvir/10 
to approximate the density profile (see long-dashed lines in 
Fig.|SJ, and equation JSJ with /3 v i r = 0.8 to approximate the 
anisotropy profile /3(r). The radial velocity dispersion calcu- 
lated from © is converted to o"3d and plotted (filled squares) 
in Figure |S] It is very close to the measured velocity disper- 
sions which confirms our expectation that high-a material is 
a more concentrated and colder subset of particles in dynam- 
ical equilibrium within the total dark matter potential. This 
is just the opposite situation to the one of surviving CDM 
subhaloes, which are a more extended, hot subset in equi- 
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Figure 6. Distribution of peaks above 2.5a in present-day galaxy-size haloes. Line styles are as in Fig. [3] The distribution resembles 
that of the Milky Way stellar halo. In the solar neighborhood (r ~ 0.03r v i r ), the 3D velocity dispersion is of order the circular velocity 
[V c (r) /V c ,msx. plotted as dotted line in the lower left panel] and there is a clear radial anisotropy, with j3 ~ 0.4 (see also Fig. 1101 . 



librium with the host potential ijDiemand. Moore fc S tadcl 
2004b). We also calculated an isotropic [13 (r) = 0] and a 
constant (3(r) = 0.5 model (see van den Bosch et al. 2004 
for the corresponding solution of eq. |1J but both cases do 
not to reproduce the measured 030 (r) very well (see circles 
vs. long-dashed line in the lower left panel of Fig. . 

As a consistency check we also solve eq. © for the 
entire dark matter halo, using (a,/3 v ,y) ~ (1,3,1.2) and 
/3 v ir = 0.45 to approximate the density and anisotropy pro- 
files. The resulting <T3d is plotted (open squares) in the lower 
left panel of Figure^] The agreement with the measured dis- 
persion is very good. As expected our average system from 
four isolated and relaxed galaxy haloes is very close to a 
stationary equilibrium solution. 



3.2.2 A fitting function for f3(r, v) 

Equation JSJ is not a very good approximation to the aver- 
age anisotropy (r) of the four G0-G3 p arent galaxies. The 
function used bv lMamon fc Lokasl i2005l) 



(3{r)=p a ^— (7) 
r + r a 

fits the data quite well for j3 a = 0.45 and r a = 0.065. For 
the high-a particles the radial anisotropy is larger and sets 
in further inside. This behavior is approximated by simply 
using the scale radii of the high-a density profiles © instead 
of r a — 0.065, and an amplitude /3 a that grows with v. 

f3(r,v) = j3au — - — , r av = r s /f u , 
r + r av 

/„ = exp(iV2) , p a » = 1 ~ 0Au-° 5 . (8) 

These anisotropy profiles are plotted with open triangles in 
the lower right panels of Figures 2] to [5] They also approx- 
imate the inner /3(r) of the haloes D12 and G4. Outside 
10% of the virial radius both the cluster D12 and the small 
galaxy G4 deviate substantially from the average /3(r) found 
for G0-G3. 
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Table 3. Density profile parameters of 4, 3.5, 3, 2.5, 2, and 
lcr fluctuations in parent haloes at z = 0. The profiles are ap- 
proximated with the general function in eq. Q using (a, j3, 7) = 
1.2) and a concentration c„ which is higher than the halo 
concentration c by some factor /„, i.e. c„ = /„c. The outer profile 
slopes /3„ and concentration factors /„ given here are calculated 
from the empirical formulae in These profiles are plotted in 
the corresponding figures as triangles. R e is the effective radius 
of the deprojected R 1 / 4 model (eq. [3) and r\i% is the average half 
mass radius of the high-cr material at 2 = 0. Both R e and ri« are 
given in units of the scale radius of the entire dark matter halo 
(r v j r /c). Values in brackets are for the cluster halo D12 which has 
a concentration of c = 4.5. For comparison, the bottom line gives 
the same parameters for the entire halo. 



V 


fu 




Re 


ri/2 


4 


7.39 


5.39 


(0.067) 


(0.1) 


3.5 


5.76 


4.93 


0.077 


0.12(0.18) 


3 


4.48 


4.51 


0.14 


0.28(0.28) 


2.5 


3.49 


4.13 


0.32 


0.51 


2 


2.72 


3.79 


0.71 


0.70(0.83) 


1 


1.65 


3.26 


2.9 


1.8(1.2) 


halo 


1.0 


3.0 


8.3 


3.5(1.8) 



3.2.3 Density slope-velocity anisotropy relation 

lHansen fc Moorel <l2004l) have recently found a relation be- 
tween the logarithmic slope of a density profile and the its 
velocity anisotropy: /3 = 771 — 772 (dlogp/dlogr), with rji 
ranging from -0.45 to 0.05 and rj2 from 0.1 to 0.35. This 
relation approximates values measured in a variety of equi- 
librium N-body and SPH systems, including CDM haloes. 
It is interesting to ask whether a similar relation also ex- 
ist for our high-cr subsets of present-day galaxy haloes. By 
combining the fitting functions for p(r, v) (eq.|2J and (3(r, v) 
(eq. |5J one indeed finds the same simple relation, 



0(d log pi d log r) = — rj 2 (7 + d log p/ d log r) 



(9) 



where 7 = 1.2 is the slope of the inner density cusp (p oc 
r -1 ' 2 ). The coefficient 772 depends only mildly on v. Parent 
CDM haloes and low-cr subsets have 772 — 0.3, while high-cr 
subsets have smaller values, r/2 = (0.29, 0.28, 0.23, 0.19) for 
v=(l, 2, 3, 4). The scatter around these relations is about 
0.1 when /3 < 0.5 (i.e. in the inner regions), and is larger 
(up to 0.4) in the outer parts. A larger, more representative 
sample of CDM haloes should be studied to quantify more 
accurately the exact slope and scatter of this important re- 
lation. The fact that even the high-cr subsets follow a nearly 
universal relation between velocity anisotropy and density 
slope supports the existen ce of a fundamental co nnection 
between the two quantities iHansen fc Moorel2004l) . This re- 
lation breaks the mass-anisotropy degeneracy present when 
one uses the line-of-sight velocities of extragalactic stellar 
halo objects (such as globular clusters) to infer total dy- 
namical masses. 




Figure 7. Shape profiles of high-cr subsets at 2 = averaged 
over four galaxy haloes. The shapes are plotted for 3.5<r (groups 
with virial mass ^ 4.9 X 10 7 Mq at z = 13.7, squares), 3o~ 
5.9 X 10 6 M Q at z = 13.7 (long-dashed lines) and > 2.2 X 10 s Mq 
at 2 = 10.5 (short-dashed lines)] and 2cr material (^ 4.8 X 1O 7 M0 
at 2 = 7.5, dotted lines). The axis ratios of the entire dark matter 
halo are plotted with solid lines. 



3.3 Shape of high-cr material 

The axis ratios a/c and b/c ar e calculated with TIPSY 1 
using the technique described in lKatd {l991). First the iner- 
tia tensors of the particles within a sphere of a given radius 
is calculated and diagonalized, then the same is done iter- 
atively for particles within the triaxial shape found in the 
previous step, until the procedure converges. Note that the 
shape at a given radius depends on all particles within this 
radius, but for density profiles less steep than p(r) oc r~ 4 
the outer particles dominate over the inner ones in the con- 
tribution to the inertia tensor. 

The average axis ratios of the four galaxy haloes G0-G3 
at the virial radius (a/c ~ 0.6, a/b ~ 0. 7) are close to the 
mean values of large samples of haloes (fj ing fc S utdl2002l : 
iFaltenbacher et aljBooS iBailin fc Steinmet j|2005l) . Figure 
!3.3l shows the average shapes of 2, 3, and 3.5cr material at the 
present epoch. Both axis ratios become smaller for higher-cr 
material at all radii. (This is also true for each of the four 
galaxies individually, not just for the average.) This result 
is consistent with the finding of Ijing fc Sutol J2002T) that 
the average axis ratios decrease with redshift, although our 
2a haloes formed at highe r redshift an d lowe r virial mass 
than the range probed bv Ijing fc Sutol i2002tl . The origin 
of the extremely prolate shape of the high-cr subsets (a : 
b : c ~ 3 : 1 : 1 for the 3.5cr material) could be a series 
of correlated head-on (i.e. low-angular momentum) mergers 
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Table 4. Mean fractional mass in present-day parent haloes of different size (2 X 10 11 Mq, 10 12 Mq, 10 13 Mq, and 10 14 Mq) contributed 
by a fixed progenitor minimum mass/rcdshift (M m ; n ,2;) selection. Only parents hosting at least one selected progenitor are included when 
computing average values and scatter. The number of parent haloes (out of ten) with selected progenitors is given in square brackets. 



Af min ,Z 


10 10 Mq, 7.0 


10 10 M Q , 4.3 


10 10 Mq, 3.1 


10 11 M , 4.3 


10 11 Mq, 3.1 


10 11 Mq, 0.8 


V 


3.2 


2.1 


1.6 


2.75 


2.07 


1.0 


2 X 10 11 Mq 


0.0 [0] 


0.049 [3] 


0.158±0.075 [5] 


0.0 [0] 


0.0 [0] 


0.567±0.071[6] 


10 12 M q 


0.0 [0] 


0.047±0.035 [8] 


0.149±0.068 [10] 


0.0 [0] 


0.123 [3] 


0.598±0.074[10] 


10 13 M o 


0.0033±0.0028 [8] 


0.077±0.028 [10] 


0.183±0.040[10] 


0.036±0.025 [9] 


0.111±0.045 [10] 


0.539±0.077 [10] 


10 14 Mq 


0.0037± 0.0021 [10] 


0.072±0.014 [10] 


0.179±0.022 [10] 


0.032±0.011 [10] 


0.116±0.023 [10] 


0.549±0.053 [10] 



along a filament aligned with the long axis jMoore et alJ 

For the total dark matter halo the axis orient ations 
change with radius, as found in IJing fc Sutol (|2002f> ■ and 
a similar behavior is observed for the high-a subsets. The 
highest a subsets show a slightly better alignment which is 
simply explained with the fact that inner particles dominate 
the shape calculation due to the very high concentrations of 
these subsets. Between the different subsets the axes are 
generally well aligned only when they have many particles 
in common, i.e. in the inner parts. As the fraction of high-a 
material drops with increasing radius, the alignment with 
the shape of the whole halo becomes worse. 



3.4 Mass fraction of high-a material 

It is well known that high-a peaks at early times are bi- 
ased towards overdense regions where larger haloes will form 
later (see Figs. and But how large is the mass fraction 
of such peaks within present-day haloes, and how does this 
fraction depend on the mass of the z — parent host? Do 
haloes which correspond to higher-a peaks today (like mas- 
sive clusters) contain a larger fraction of early high-a ma- 
terial? This question is not well defined, since the fraction 
of particles from peaks above 2a in a present-day galaxy 
halo grows when the peaks are selected with a larger mass 
threshold at lower redshift (see upper right panel of Fig. 
0. Here we use fixed mass threshold/redshift pairs to select 
progenitors and compare the mass fractions they contribute 
to parent haloes of different sizes. The selection of a fixed 
progenitor mass may be motivated in studies of old stel- 
lar populations if, for example: a) Population III stars form 
in "minihaloes" above a molecular cooling mass, i.e. massive 
enough to allow efficient gas cooling via roto- vibrational lev- 
els of H 2 , M > Mh 2 « 6 x 10 5 [(1 + z)/20]" 3/2 M Q (virial 
temperatures > 2000 K); and b) Metal-poor halo stars (Pop- 
ulation II) and globular clusters form in haloes above an 
atomic cooling mass, i.e. massive enough to allow efficient 
gas cooling and fragmentation via excitation of hydrogen 
Lya, M > Mh » 10 8 [(1 + 2)/10]~ 3/2 M Q (virial tempera- 
ture > 10 4 K). 

To test whether, at the present epoch, massive clus- 
ters contain a different fraction of high-a material than field 
galaxies, we have used a lower-resolution N-body simulation 
(300 3 particles in a 90 Mpc box, for a resolution of 10 9 Mq) 
in order to obtain a larger galaxy /cluster sample. We have 
selected ten parent hosts at z = with FOF masses close 
to 2 x 10 11 M©, 10 12 M , 10 13 M Q , and 10 14 M , for a to- 
tal of fourty haloes. We have then marked FOF groups with 



more than 100 particles at redshifts 7.0, 4.3, 3.1 and 0.8, and 
determined the fractions of the virial mass of the parent be- 
longing to a fixed progenitor mass/redshift (A/, z) pair. The 
values obtained are given in Table 01 the mass fractions are 
practically constant for all parents which host any of the 
selected progenitors. As we select parents of lower masses, 
the number of parents hosting rarer (higher-a) progenitors 
drops from 10 to within about a decade in mass. 



4 SOME APPLICATIONS 

We have shown that the final distribution and kinematics 
of dark matter particles selected from early branches of the 
merger tree are systematically different than those of the 
parent halo as a whole. These properties are also relevant 
for old stellar populations if these form predominantly in 
early low-mass progenitor haloes, as stars behave essentially 
as collisionless systems just like the dark matter particles in 
our simulations. In the following we briefly discuss a number 
of possible applications. 

4.1 Remnants of the first stars 

Numerical simulations performed in the context of hierarchi- 
cal structure formation theories suggest that the first (Pop- 
ulation III) stars may have formed out of metal-free gas in 
dark matter minihaloes of mass above 6 x 10 5 M (Abel et 
al. 2000; Bromm et al. 2002; Yoshida et al. 2003; Kuhlen 
& Madau 2005) condensing from rare high-a peaks of the 
primordial density fluctuation field at z > 20, and were 
likely very massive. Barring any fine tuning of the initial 
mass function (IMF) of Population III stars, intermediate- 
mass black holes (IMBHs) - with masses above the 5-20 M 
range of known 'stellar-mass' holes - may be one of the in- 
evitable endproduct of the first episodes of pregalactic star 
formation (Madau & Rees 2001). 

Where do relic pregalactic IMBHs lurk in present-day 
galaxy halos? To shed some light on this question, we have 
populated our 3a (3.5a) simulated progenitors at zo.9 = 17.9 
(2:0.9 = 21.2) with one seed IMBH for every 6 x 10 5 so- 
lar masses of halo material. As discussed by Volonteri et al. 
(2003), these IMBHs will undergo a variety of processes dur- 
ing the hierarchical buildup of larger and larger haloes, like 
gas accretion, binary hardening, black hole mergers, triple 
interactions. While we neglect all of these effects here, our 
dark matter simulations do correctly model the bias in the 
formation sites, the accretion into larger haloes, and the 
competing effects of dynamical friction and tidal striping 
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Figure 8. Distribution of peaks above 3<j in present-day galaxy-size haloes. Line styles are as in Fig. [3] The circles in the upper left 
panel show the estimated mass and number density of wandering IMBHs assuming J\T = 1 seed hole for every 6 X 10 5 solar masses of 3<r 
progenitor halo material at 20.9 = 17.9, and a final black hole mass of 1.5 X 10 4 Mq (see text for details). 



within larger potential wells, as these complicated dynam- 
ical processes are dominated by the dark haloes that host 
the black holes. Therefore, in our toy model, the distribu- 
tion of 3 or 3.5a material at z = describes the properties 
of holes wandering through within today's galaxy haloes. 
The predicted IMBH number density and mass density pro- 
files are shown as circles in Figures |H] and |5] the former may 
be regarded as an upper limit since we have neglected black 
hole mergers, while the latter have been estimated assuming 
that these off-rmiclear black holes have grown by accretion 
to a mean mass of 1.5 x 10 4 Mfr>, which is a rough estimate 
obtained from Figure 14 of IVolonteri et all i200St) . 

Depending on the IMF of Population III objects, some 
first-generation low-mass stars may have survived until to- 
day. Their number density profile n(r) within the Milky Way 
can again be read following the circles in Figure |H1 under 
the assumption that N = 1 metal-free star survives for ev- 
ery 6 x 10 s solar masses of 3cr progenitor halo material at 
zo. 9 = 17.9. It is easy then to scale up the predicted value 
of n(r) if N 3> 1 such stars were to survive instead. On 
average, we find that about 1/3 of these remnants would 



lie today in the bulge, i.e. in the inner 3 kpc. This frac- 
tion fluctuates between 24% and 45% in the four galaxies 
G0-G3. The density in the solar neighborhood is of order 
~ 0.1 TV kpc -3 , three orders of magnitude lower than in the 
bulge. The number density of remnants would be lower and 
more concentrated toward the galactic center if Population 
III stars only formed within rarer 3.5cr peaks (Fig. |5J. On 
average, about 59% of them would now lie within the bulge 
(range is 38%-84%) and the local number density would be 
only ~ 0.02 N kpc -3 . Even for N = 10 - 100, this is an 
extremely small value, many orders of magnitude below the 
local number density of halo stars. The above results suggest 
that the very oldest stars and their remnants should be best 
searched for within the Milky Way bulge. 



4.2 Stellar haloes 

Material from > 2.5cr peaks has today a density profile that 
is very similar to t he stellar halo around the Milky Way 
jMoore et al. 1 ^005 ). It contributes a few percent of the to- 
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Figure 9. Same as Fig. HI but for peaks above 3.5cr and for Population III progenitors collapsing at 20.9 = 21.2. 



tal virial mass, and therefore contains enough baryons to 
build up a 10 9 Mq stellar halo with a reasonable star for- 
mation efficiency. The assumption of a common pregalac- 
tic origin between such a stellar component and the sur- 
viving Local Group dwarf galaxies provides an additional 
constraint and allows us to determine the progenitor mass 
threshold/red shift pai r whic h best fit the data. From this ar- 
gument iMoore et al~l (|20nfJ) identified hosts above 10 8 M Q 
at z — 12 as the progenitor haloes which the bulk of halo 
stars originally belonged to. 

To check whether this simple model reproduces the 
kinematics as well as the radial distribution of halo stars, we 
have compared the predict ed radial velocit y dispe rsion pro- 
file with recent data from iBattaglia et alJ (1200.4) (see Fig. 
I10II . All haloes were rescaled (i.e. r v i r — ► fr v - a , V c — > fV c , 
M v i T — > / 3 M v j r ) to produce a local circular velocity of 
220 km/s after taking into account the increase in cir- 
cular velocity due to bar yonic contrac t ion [1 4(8.5 kpc) = 
1.125 V c ,max,DM] found bv lMaccio et alJ i2005t) for our halo 
Gl. The rescaled virial masses are given in Figure [TU1 The 
predicted er r (r) are close to the observed ones, especially for 
halo GO and G4. Our model predicts a radial anisotropy that 



grows with radius (see lower right panel of Fig. 1101 : in the 
solar neighborhood it gives 03D — V e and ~ 0.4, in goo d 
agreement with the observed values dChiba fc Beersl F2000). 
The declining velocity dispersion near the center is charac- 
teristic of haloes containing only dark matter. In the Galaxy 
we expect a flatter profile in the inner, baryon-dominated re- 
gions. 

Both a high concentration (c = 18), highly anisotropic 
halo like G4 and a c = 10 halo with smaller anisotropy like 
GO fit the data equally well. A detailed comparison requires 
knowledge of the tangential components of the velocity dis- 
persion, in order to determine the concentration and virial 
mass of the Milky Way halo. The recently discovered stellar 
halo of M31 has a similar density profile than that of our 
Galaxy iGuhathakurta et al.l l2005'l . The line-of-sight veloc- 
ity dispersion observed in M31 is mostly due to tangential 
motions: preliminary results yiel d a constant at I y2 of abou t 
80 km/s for 50 < r < 150 kpc iGuhathakurta et alJkoOEfi . 
close to the value expected for the GO scenario but much 
higher than that of the G4 model. More data for M31 will 
hopefully soon become available, and many radial velocities 
of Milky Way halo stars will be obtained by surveys like 
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Figure 10. Kinematics of our model stellar haloes compared with observations of Milky Way halo stars. The simulated CDM haloes are 
rescaled to give a local circular velocity of 220 km / s after taking baryonic contraction into account (see text for details) . The lower right 
panel shows the predicted anisotropy profiles fi(r) = 1 — o- 2 (r)/2cr 2 (r) for the G0-G3 average halo (dashed line: halo 's tars', solid line: 
dark matter) and for halo stars in GO (dash- dotted) and G4 (dotted). The data point is the local halo anisotropy from lGniba fc Beersl 
i200(l) . The curves in the other three panels show, from bottom to top, the 'stellar' at / y/(2) ~ ag ~ cr^ (dotted line), oy (solid line) and 
(dash- dotted line). The uppermost solid curve depicts the 3D velocity dispersion profile of the dark matter (note that it is much larger 
and smoother than the stellar <T3d). To illustrate the fluctuations in the stellar dispersions, ay and 030, dm in halo GO are plotted twice: 
once using our usual choice of 30 logarithmic bins to the virial radius, and once using smaller linear bins (open circles with Poissonia n 
error bars for ay, another solid line for <T3d,dm)- The points are the measured oy of Milky Way halo stars from iBattaglia et alJ (2005), 
with (squares) and without (stars) the contribution from satellite galaxies. The error bars reflect only the accuracy of the observations, 
while Poisson noise (which causes similar, in some bins even larger, errors) is not included. The declining dispersions near the center are 
characteristic for halo containing only dark matter. In the Galaxy we expect a flatter profile in the inner, baryon-dominated regions. 



RAVE and SDSS-II. In the next decade the planned GAIA 
satellite will provide very detailed phase-space positions. 

Another interesting feature of our model is that it pre- 
dicts significant deviations from the simple radial velocity 
dispersion profiles expected for a smooth stellar halo i n an 
NFW potential (like the models in lBattaelia et, all200Bl. Fig. 
4). These deviations are evident in the upper right panel of 
Figure 1101 (but not in the lower left panel where we show 
an average over four haloes). 2 The fluctuations are not an 

2 The sample from Battag fia et alJ l2005h and our G4 outer stel- 
lar halo contain similar numbers of 'stars'. Both have relatively 



artifact of our analysis: with our usual choice of 30 loga- 
rithmic spherical bins out to the virial radius [similar to 
the binning of the observed oy (r)] the fluctuations often ex- 
tend over various bins with each bin containing thousands 
of 'stars'. To show fluctuations of oy(r) on smaller scales we 
have binned 'stars' in halo GO using linear spherical bins, 
and have ploted the resulting oy(r) profile in the same ra- 

large Poisson noise and are consistent with smooth oy(r) profiles. 
In haloes resolved with more particles like G0-G3 and D12, the 
'stellar' radial velocity dispersions are clearly inconsistent with 
smooth profiles. 
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dial range of the observations (circles with error bars in the 
upper right panel of Fig. [TUl . The linear bins contain at least 
108 stars, the error bars represent Poisson noise. We have 
also plotted o"3d,dm(V) of halo GO twice using different bin 
sizes. Note that the velocity dispersion of the dark matter 
halo as a whole (instead of the material belonging to the 
> 2.5a z = 12 subset) in our simulated galaxies is generally 
a smooth function of radius, with small deviations becoming 
apparent only at r ~ 200 kpc. As fluctuations in the stel- 
lar oy(r) are not caused by mass bound to subhaloes (the 
mass fraction is much too small particularly within the inner 
halo), this leaves tidal streams as the most probable cause 
of such deviations. The entire dark halo is assembled from 
about 10 5 resolved progenitors, the tidal debris of which are 
well mixed and smooth. In our model, the stellar halo is 
made up instead from only a few hundred building blocks 
(and a significant fraction of it comes from the most massive 
ones), and is therefore expected to show some granularity. 
A more detailed study of tidal features in our model stellar 
haloes is left to a future paper. 

4-2.1 Comparison to recent related models for the 
formation of the stellar halo 

During the completion of this work two related models 
for the formation of th e stellar halo iBullock fc Johnston! 
120051: lAbadi et al]l2005fi came to our attention. Both stud- 
ies are more detailed than our in the sense that they try 
to mo del the effect of star form ation and chemical evo- 
lution. iBullock fc Johnston! l|2Q05) use a semi-analytic ap- 
proach combined with N-body simulations of satellite dis- 
ruption in an external, growing galaxy potential. Like in this 
work the number of luminous satellites and the mass of the 
final stellar halo is adjusted to match the observations by as- 
suming that reionisation suppresses star formation in later- 
forming small progenitors. Our stellar haloes, however, seem 
to be more concentrated than those of IBullock fc Johnston! 
( 2005). This may due to the fact that we only tra ce old stel- 
lar populations, while IBullock fc Johnstonl i2005T) also allow 
for more recent episodes of star formation in their proto- 
galactic building blocks. They also model the disruption of 
satellites assuming an external spherical potential: this is 
probably a good approximation to recent accretion events 
which build up the outer stellar halo. The inner stellar halo, 
however, is built up through a series of early, major merg- 
ers: these cause rapid potential fluctuations (violent relax- 
ation) which are not accounted for in the spherical poten- 
tial approximation. This li mitation was also pointed out by 
IBullock fc Johnstonl <2005ft . 

The cosmological SPH simulations of lAbadi et alJ 

(2005), while not including any global mechanism for sup- 
pressing star formation in later-forming small progenitors, 
do lead to similar stellar halo profiles and kinematics as our 
model. They seem to produce, however, far too many halo 
stars: recent SPH simulations at higher resolution and with- 
out strong feedback effect s also overpredict th e number of 
satellites by a large factor iMaccio et alJ <l2005f) and tend to 
produce even more halo stars. The similarity of the pro- 
file shape and kinematics of stellar haloes in all these mod- 
els suggests that they are a robust and g eneric outcome of 
hierarchical structure f ormation (see also lHansen fc Moorel 
120041: iDekel et afll2005l) and do not depend on the detailed 



formation history. Note that this does not contradict our ar- 
gument about the possibility to use the radial extent of col- 
lisionless tracer populations to learn about their formation 
epoch and sites. It just illustrates how this arguments is lim- 
ited to fossil records from very high redshift (z > 8, see also 
§ 12.21 . At lower redshifts the typical haloes are larger and 
the exact star formation sites must be taken into account. 
Efficient cooling (and/or numerical losses of gas angular mo- 
mentum to the dark matter in under-resolved progenitors) 
can produce a population of stars at some intermediate red- 
shift with a distribution that is similar to much older pop- 
ulations from early high-cr peaks. By z — the two sub- 
sets would have very similar concentrations and kinematics, 
which might expla in the similar results found here and in 
lAbadi et all J2005D . 

4.3 Metal-poor globular clusters 

The clustering properties of metal-poor globular clusters 
contain clues on their formation sites and the epoch when 
star formation was suppressed by feedback processes (e.g. 
reionisation, supernova-driven winds) in low-mass haloes. In 
the Milky Way metal-poor globular clusters follow the same 
radial profile as halo stars, suggesting within the framework 
of our model a common origin within early 2.5cr peak pro- 
genitors at z « 12. Observations of the distribution of glob- 
ular clusters within present-day haloes of different masses 
could provide infor mation on feedback effects as a func- 
tion of environment jMoore et al. 1 120051) . The suppression 
of globular cluster formation at some early epoch may also 
explain the bimoda lity observed in cluster metallicities (e.g. 
IStrader et aTll2005h . The widely used assumption that glob- 
ular cluster formation is a fair tracer of star formation, com- 
bined with the suppression of the formation of metal-poor 
globulars after reionisation, imply that the amount of high-cr 
material in a halo is proportional to the number of metal- 
poor globular clusters. From the results in § 13.41 it follows 
then that a simple universal reionisation epoch would lead 
to a constant abundance of metal-poor globular clusters per 
virial mass. Deviations from this simplest case may provide 
information about the local reionisation epoch, and whether 
regions with more (less) metal-po or globulars per viri al mass 
were reionised later (earlier) fsee lMoore et al~ll2005l) . 

4.4 Elliptical galaxies 

The projected luminosity profiles of elliptical galaxies fol- 
low the R 1 / 4 law and resemble re scaled versions of pro- 
jected CDM halo density profiles jLokas fc Mamonl |200ll : 
iMerrit et al1l2005ft . similar to the 2 = profiles of our high- 
er subsets. Ordinary elliptical galaxies may form by multiple 
mergers or by th e major merger of two disk galaxies (e.g. 
Dekeletal . 2005). Early formation (z ~ 6) in a large num- 
ber of progenitors followed by a series of gas-poor, essentially 
collisionless mergers could build up the giant ellipticals (cD) 
in the centers of galaxy clusters. These mergers might undo 
some of the early dissipationa l contraction of the total mass 
distribution <Gao et al]l2004Tl . 

Figure fTTI shows the distribution of material belonging 
to different a peaks in a z = cluster halo. The inner cluster 
regions are entirely made up of material from > 2a peaks, 
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Figure 11. Distribution of peaks above 1 (short-dashed), 2 (long- dashed), 3 (dotted) and 4a (dash-dotted) in a present-day cluster-size 
halo (M vir = 3.1 X 10 14 M ). The redshift and minimum mass pairs used are (z = 1.6, 7.4 X 10 9 Mq) for the la peaks, (4.4, 6.3 X 10 9 Mq) 
for the 2a, (7.4, 4.3 X 10 9 Mq) for the 3a, and (8.7, 1.8 X 10 10 Mq) for the 4a peaks. The triangles in the upper left panel show the 
fitting function l2*t . 



and this selection contributes about 10% to the cluster virial 
mass. The central brightest galaxies in SDSS clusters have 
an average effective radius of R e ~ 20 kpc and presumabl y 
about 10% of the clusters luminosity (Zib etti et al.l 12005). 
Our 2a subsets contain enough mass to account for the ob- 
served cD luminosities but they are more extended. If the cD 
stars formed early in the inner parts of all progenitor haloes 
above 2a then they would be found closer to the cluster cen- 
ter today (see Section f2.2H . A simple way to select the inner 
parts of 2a progenitors is to mark higher-cr peaks at earlier 
epochs. Our 3a selections have a realistic effective radius of 
about 20 kpc. Therefore we predict the kinematics of the 
average brightest cluster galaxies to follow the 3a selections 
shown in Figure 11 1 1 1 if (or where) the real cluster potentials 
are similar to the uncontracted potentials formed in dissipa- 
tionle ss cosmological simulations (as suggested bv lGao et alJ 
12004) . 

The anisotropy of a z — subset depends mostly on 
its density profile shape (or v) and not on the selection red- 
shift or the number of mergers it has undergone (Figures 



2Jto[5Jl A few mergers are enough to establish the general 
(3(r) profile JSJ and the slo pe - anisotropy rela tion (Sec- 
tion 13.2.31 . This agrees with Mo ore et alJ B2004i) who find 
(3 ~ 0.5 after only one low angular mo mentum merger of 
initia l ly isotropic CDM haloes (see also lHansen fc Moorel 
12004) . foekel et alJ J2005I) also find a similar f3(r) profiles af- 
ter only one merger, which even starts with disk galaxies, 
i.e. nearly circular initial orbits for the sta rs (see their (3 
for the old stars, Fig. 1 of lDekel et aDl2005T) . Their average 
anisotropy at the half mass radius of the old stellar compo- 
nent is P = 0.3, which is lower than our average of /3 = 0.45, 
but considering the very different initial conditions of the 
stars and the larger numbers of mergers in our cosmological 
context the results are surprisingly similar. Therefore we ar- 
gue that our collisionless results are relevant both for giant 
elliptical gala xies which form ed in a series of dissipation- 
less mergers (iGao et alJl2004) and also for the old stars in 
smaller ellipticals. 

The general relation between density pofile slope and 
anisotropy for high sigma subsets /3(dlogp/cilogr) ~ 
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Figure 12. Same as Fig. |l]for peaks above 1 (short- dashed), 2 (long- dashed), 3 (dotted) and 3.5c (dash-dotted) in a small galaxy halo 
(M vir = 10 11 M s ) at z = 0. The redshift and minimum mass pairs used are (z = 2.4, 4.8 X 10 s Mq) for the la peaks, (8.7, 7.6 X 10 6 Mq) 
for the 2a, (10.5, 2.2 X 10 s Mq) for the 3<r, and (13.7, 4.9 X 1O 7 M ) for the 3.5<r peaks. The triangles in the upper left panel show the 
fitting function l2*t . 



— 0.23(1. 2 + d log p/dlog r) can be applied to elliptical galax- 
ies to infer f3(r) from the observed stellar profile. This is use- 
ful to obtain dynamical mass estimates for ellip tical galax- 
ies, w here one has to assume some /3(r) iMamon fc Lokasl 
2005) and it is probably the best guess before realistic, cos- 
mological, hydrodynamical simulations of elliptical galaxy 
formation become available. 



5 CONCLUSIONS 

We have used high-resolution cosmological N-body simu- 
lations to trace the spatial distribution and kinematics in 
present-day CDM haloes of collisionless material that origi- 
nally belonged to selected branches of the merger hierarchy. 
Our main results can be summarized as follows: 

(i) Hierarchical merging does not efficiently mix haloes 
up. The final distribution of dark matter particles retains a 
memory of when and where they collapsed initially, allowing 



a unique test of the popular bottom-up paradigm for the 
formation of cosmic structure. 

(ii) Today's distribution and kinematics of halo substruc- 
ture that formed at very high redshift depends mostly on 
the rarity of the primordial density fluctuations they corre- 
spond to. For example, the z — density, shape, anisotropy, 
and velocity dispersion profiles of material originating from 
early 3a peaks is independent of the redshift/minimum mass 
pair used to select it. The mass fractions within the parent 
host, however, grow if such particles are choosen from lower 
redshift progenitors at a higher mass threshold. 

(iii) High-cr material should be looked for close to the cen- 
ter of massive parent haloes. The concentration and outer 
slope of the mass density profile are larger for the rarer 
peaks. 

(iv) The anisotropy increases faster with radius for ma- 
terial originating within rarer peaks (eq.|HJl. The parameter 
f3(f) steepens with increasing v at the same rate as the half- 



mass radius shrinks, i.e. we generally find f3(r 



1/2) 



0.45. 



(v) Velocity dispersions are lower for material from the 
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rarer peaks, and particle orbits are more radial. While the 
average velocity dispersion profiles agree with stationary so- 
lutions to the Jeans equation (gj, the profiles of individual 
haloes show significant structure and are inconsistent with 
smooth a(r) solutions (see Fig. If 01 . 

(vi) The high-er subsets have much more elongated shapes 
than their host haloes. For 3a material the mean c/a is about 
0.35, while present-day galaxy-size haloes have a mean of 
0.6. 

(vii) These properties do not depend on which regions 
within early (z > 8) progenitors are marked and traced to 
the present time. Particles originating within the cores of 
high-CT progenitors are distributed at z — in the same way 
as the entire virial mass of these early structures. This is not 
true at later epochs, where the marking of only the densest 
progenitor regions results in more concentrated density pro- 
files (e.g. the selection of only the densest f 0% of the mass 
of progenitor halos above la at z — 1.6 yields a present-day 
distribution similar to that of all particles above 2a). 

(viii) If the first stars form at early epochs in peaks above 
3.5a, then half of their remnants should be found in the 
bulge, within 3 kpc of the galactic center. In the solar neigh- 
borhood the density of such very old population is 1000 
times lower than in the bulge. Also, their characteristic ve- 
locities are 2.5 times lower than those of dark matter parti- 
cles, and their orbits more radial (/? ~ 0.5). 

(fx) The radial profile of the stellar halo and metal-poor 
globular clusters of the Milky Way suggest that these com- 
ponents formed in rare early peaks above 2.5a at redshift 
above 10. Typical outer halo objects have radial orbits, and 
become isotropic near the galactic center. 

(x) Radial orbits are a general outcome for any concen- 
trated stellar component assembled through gas-poor merg- 
ers. The anisotropy parameter (3(r,v) correlates well with 
p(r, v) and is not sensitive to the detailed assembly history 
of such a component. 

The applications discussed above should be taken as a 
first attempt at predicting the spatial distribution and kine- 
matics in present-day galaxies of objects that formed within 
early protogalactic systems, all in the context of hierarchical 
structure formation theories. Some interesting directions for 
future work may include, e.g., combining spatial information 
on Population III remnants with semi-analytic prescriptions 
for the growth and dynamical evolution of black holes, as 
well as using chemical evolution models to translate the age 
gradient for stellar haloes found here into metal abundance 
gradients. The techniques presented here also allow to study 
the properties of stellar streams in realistic, triaxial, clumpy 
CDM galaxy haloes. 
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